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August  8,  1988 

1.  Introduction 

T  his  note  concerns  new  eigenfunction  expansions  for  the  fields  in  an 
irregular  acoustic  waveguide.  The  expansions  are  associated  with  a  non- 
selfadjoint  eigenvalue  problem  that  has  the  eigenvalue  in  some  of  the 
boundary  conditions.  The  waveguide  consists  of  two  layers  of  fluid  that 
are  separated  by  a  frictionless  free  interface.  The  field  expansions  con¬ 
verge  uniformly  on  both  sides  of  the  interface  even  if  the  fluid  density  is 
discontinuous  at  the  interface. 

We  use  the  eigenfunction  expansions  to  convert  the  reduced  wave 
equation  and  the  conditions  at  the  boundaries  of  the  waveguide  into 
a  system  of  first-order,  ordinary  differential  equations.  These  coupling 
equations  are  similar  to  Shevchenko’s  equations  for  coupling  between 
local  modes  in  the  waveguide  (Sov.  Phys.  Acoust.  7(1962)  pp.  392-397, 
eqs.  (16-18)). 

The  eigenvalue  problem  and  the  field  expansions  are  described  in  the 
next  section.  The  coupling  equations  are  stated  in  the  last  section. 

2.  The  eigenfunction  expansions 

A.  Notation  and  assumptions.  The  discussion  in  this  note  is  lim¬ 
ited  to  two-dimensional  problems.  The  real  variables  x  and  z  denote, 
respectively,  the  coordinates  of  horizontal  and  vertical  position  in  the 
waveguide.  The  following  notation  is  used  for  the  indicated  limits  and 
jump  discontinuities  of  an  arbitrary  function  F(x,z ): 

(1)  F(x,  z  +  0)  =  lim  Fix,  z  -F  e), 

t— >o 
e>0 

(2)  F(x,  z  -  0)  =  Jim  F(x,  z  -  t), 

*>o 


jump|F)  =  F(x,  z  +  0)  -  F(x,  z  -  0). 

(«.*) 


(3) 


Surface  2  =  0  (p  =  0) 


Upward  normal  ft  oc  1] 


Interface  z  =  —H(x) 


.  1  dp  .  . 

(p  and  -  —  continuous) 
p  on 


Bottom  2  =  —L  (p  =  0) 


Fig.  1.  Geometry  of  the  waveguide  and  conditions  at  the  boundaries. 


A  representative  section  of  the  waveguide  is  sketched  in  Figure  1. 
t  he  surface  boundary  coincides  with  the  horizontal  line  z  =  0,  and 
the  bottom  boundary  coincides  with  the  horizontal  line  z  =  -L.  The 
interface  between  layers  in  the  waveguide  is  described  by  the  curve  z  = 
-H (x).  It  is  assumed  that  0  <  H{x)  <  L  for  all  x  and  that  the  function 
H(x)  has  a  continuous  second  derivative. 

The  local  sound  speed  c(x,  z)  is  continuous  in  each  layer,  and  the  local 
mass  density  p(x,  z)  is  continuously  differentiable  in  each  layer.  These 
quantities  may  be  discontinuous  at  the  interface,  but  it  is  assumed  that 
each  of  the  four  limits  c(x,  -H{x)  ±  0)  and  p{x, -H(x)  ±  0)  exists  and 
is  positive. 

Sound  in  the  waveguide  is  assumed  to  depend  harmonically  on  time 
with  angular  frequency  w.  The  acoustic  pressure  at  time  t  can  be  ex¬ 
pressed  as  the  real  part  of  the  product  p(x,  z)  exp(-iud),  where  the  com¬ 
plex  pressure  p(x,  z)  solves  a  reduced  wave  equation  in  two  dimensions: 
if  -L  <  z  <  -  H{x)  or  -H{x)  <  z  <  0  then 


(4)  0 &(*b)  s**’ !)) + 1 ■ ’> £(*b)  i**' ) 


ID' 


C(X,  2)* 


p(x,z)  =  S(x,z). 


fit 
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The  forcing  term  S(x,  z)  in  eq.(4)  represents  distributed  sources  of  sound 
inside  the  waveguide. 

The  acoustic  pressure  vanishes  at  the  external  boundaries  of  the  wave¬ 
guide,  and  it  is  continuous  at  the  interface: 


(5) 

p(x,-L)  =  0, 

(6) 

p(x,0)  =  0, 

(7) 

jump  (p)  =  0. 

(*,-«(*)) 


The  fluid  velocity  associated  with  sound  in  the  waveguide  is  propor¬ 
tional  to  the  acoustic  pressure  gradient.  The  components  of  velocity  in 
the  horizontal  and  vertical  directions  can  be  expressed,  respectively,  as 
the  real  parts  of  the  products  u{x,z)exp(-iud)  and  tc(x,  z)exp(-»u>t), 
where  the  complex  velocity  components  u(x,  z)  and  tv(x,  z)  are 


(8) 

u(x,z)  = 

0) 

w(x,z)  = 

1  d 


iup{x,z)  dx 


P(*.  *), 


1  d  .  , 


iojp{x,z)  dz 


At  the  interface  the  component  of  fluid  velocity  normal  to  the  interface  is 
continuous.  This  requirement  can  be  expressed  in  terms  of  the  complex 
velocity  components  as 


(10) 


jump  fu>  -f  H(x)  u)  =  0, 
(*-«(*)) 


where  //(x)  =  ( d/dx)H(x ). 

B.  The  eigenvalue  problem.  At  each  position  x  the  following  equa¬ 
tions  determine  countably  infinite  sets  of  eigenvalues  kn(x)  and  eigen¬ 
functions  4>n(x,z): 


(H) 

p(x' 2)  iz  2)) +  " kn{x)i)  l)  ® 0 

if  -L  <  z  <  -H(x)  or  -H{x)  <  z  <  0, 
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(12) 


<j>n(x,-L)  =  0, 


(13) 


<^B(x,0)  =  0, 


(14) 


jump 


~ikn(x)<f> n 


=  0, 


(15)  jump  [-f^-^n  +  ^(*)«*»(*)^n)]=0- 

(*,-»(*))  1  P  \OZ  )\ 

The  eigenfunctions  are  assumed  to  be  normalized  in  such  a  way  that  the 
mixed  partial  derivatives  (d2/dx  dz)4>n(x,Q)  exist  at  the  surface. 

For  every  x  all  but  finitely  many  eigenvalues  kn(x)  are  purely  imagi¬ 
nary.  The  indices  n  of  the  eigenvalues  run  through  the  integers  0,  ±1 ,  ±2, 
. . . ,  and  they  can  be  assigned  to  the  eigenvalues  in  such  a  way  that  for 
some  integer  IV(x)  the  following  statements  are  true: 

(a)  if  |n|  >  N{x)  then  kn(x)  is  imaginary  and  |fc„(x)|  =  0(M)» 

(b)  if  n  >  N{x)  then  -iJfc„(x)  is  positive  and  increases  in  value  as  n 
increases; 

(c)  if  n  <  -N(x)  then  -ikn{x)  is  negative  and  decreases  in  value  as 
n  decreases. 

C.  Acoustic  field  expansions.  For  each  index  n  define 


(16)  a„(x)  =  J  ^  (kn(x)p(z,z)-i-^p{x,z)j  -^—^4>n(x,z)dz 


-r\-rp{x,-H{z))  jump 


1  °S 

~pd~**n" 


(!7) 

d„(x)  =  2  &n(x)  [  ~  <t>n(x,  z)2dz  + 

J-L  ?(*<*) 


1 

r—  jump 

M*)  (*,-«(*)) 


The  following  expansions  of  the  complex  acoustic  fields  converge  uni¬ 
formly  in  the  intervals  -L  <  z  <  ~H(x)  and  -H(x)  <  z  <  0: 


(18) 


p(*,z)-  E 
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(19) 


U(*,2)  = 


1 

iup(x,z) 


®«(x> 


(20) 


w(x,z)  = 


.(*)  d 


iup{x,  z)  n£^ao  d„(x)  dz 


4>n(* i  Z). 


In  particular,  expansions  (19)  and  (20)  do  not  exhibit  Gibb’s  phenom¬ 
enon  at  z  =  -H(x)  when  the  velocity  components  u(x,  z)  and  u>(x,z) 
are  discontinuous  at  the  interface. 


3.  The  coupling  equations 


The  functions  a„(x)  that  appear  as  coefficients  in  the  acoustic  field 
expansions  solve  an  infinite  system  of  first-order,  ordinary  differential 
equations.  Each  equation  in  this  system  takes  the  form 

(21)  ±an{x)  =  ikm(x)an(x)~  £  Jj=^cm„(x)  -  isn(z), 

'  m=~oo  m '  ' 

(for  n  =  0,dbl,±2,. . . ) 

where  the  forcing  term  a„(x)  is  a  moment  of  the  source  distribution: 

l 

(22)  sB(x)  =  /  S(x,z)~ -  -  <{>n(x ,  z)  dz . 

J-L  P(X,  Z) 

It  is  possible  to  express  the  general  coupling  coefficient  cran(x)  in 
terms  of  several  factors  to  which  physical  meaning  can  be  attributed. 
First,  define  the  following  quantities: 

(23)  Lmn(x,  z)  =  -j^4>m(x,  z)^$„(x,  z) 

+  ^fcm(*)kn(x)  —  <f>m{x,z)<t>„{x,z), 

(24) 

0  0 

Mmn(x,z)  =  ~<{>m{x,z)~<t>H(x,z)  -  km(x)  k„(x)  0m(x,  z)  4>n{x,  z), 


(25) 

&  0 

^m«(®i  z)  —  z)~^<f>n(x,  z)  +  km(x)  kn{x)  z)  <j>n(x,  z) 

+  N{*)*(km(x)<t>m{xtz)-~<t>n(x,z) 

0  \ 

+  fo<fim(x,z)kn(x)4n(xtz)  ), 
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(26)  U«h))-6M£  Gdb)- 

Use  these  quantities  to  construct  the  following  factors: 

(27)  Am.(.)  =  £ 

-L&irtb))  L-'(x-‘>d‘- 

(28)  Imn(x)  =  ^T(x)  jump 

(29)  J  mn  (x)  =  H(x )  jump  [Mmn/p] , 

<*,-»<*)) 


(30)  tfmB(x) 


; - ^r-rr  jump 

i  +  H(xy  (.,-«(*» 


|Afmn/p) 


+ 


ff(*)i 

- : -  jump 

1  +  H(x)2  (*,-«(*)) 


[jVmn  Dfl|l/p)]  , 


where  H(x)  =  (d2/dx2)H(x)  and  m,  n  ss  0,  ±1,  ±2, .... 

The  factor  /4mB(x)  is  associated  with  lateral  variations  in  material 
properties  of  the  fluids  inside  the  waveguide.  The  factors  Imn(x),  Jmn(x) 
and  Kmn(x)  are  associated  with  conditions  at  the  interface.  These  fac¬ 
tors  can  differ  from  zero  only  if  material  properties  of  the  fluids  are 
discontinuous  at  the  interface,  and  if  the  slope  or  the  curvature  of  the 
interface  does  not  vanish.  The  following  theorem  states  the  relation  of 
these  factors  to  the  general  coupling  coefficient  em„(x). 


Theorem.  If  m  n  then 

(31)  c„„(l)  =  ('>»•(■')  +  '—W  +  5^)  (I) 


However ,  if  m  =  n  then 

(32)  CBB(X)  =  --<*»(*)  +  Jnn(x)  +  2km(x)2  Knn^' 
where  dn(x)  =  (d/dx)dn(x). 
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